Inhomogeneity in the Supernova Remnant Distribution as the 
Origin of the PAMELA Anomaly 



Nir J. Shaviv 1 , Ehud Nakar 2 & Tsvi Piran 1 

1. Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 91904, Israel 
2. The Raymond and Beverly Sackler School of Physics & Astronomy, Tel-Aviv University, 

Tel-Aviv 69978, Israel 

ABSTRACT 



Recent measurements of the po sitron /electron ratio in the cosmic ray (CR) 
flux exhibits an apparent anomaly ( Adriani et al. II2OO9I ). whereby this ratio in- 
creases between 10 and 100 GeV. We show that inhomogeneity of CR sources on a 
scale of order a kpc, can naturally explain this anomaly. If the nearest major CR 
source is about a kpc away, then low energy electrons (~ 1 GeV) can easily reach 
us. At higher energies (> 10 GeV), the source electrons cool via synchrotron 
and inverse- Compton before reaching Earth. Pairs formed in the local vicinity 
through the proton/ISM interactions can reach Earth also at high energies, thus 
increasing the positron/electron ratio. A natural origin of source inhomogeneity 
is the strong concentration of supernovae in the galactic spiral arms. Assuming 
supernova remnants (SNRs) as the sole primary source of CRs, and taking into 
account their concentration near the galactic spiral arms, we consis tently recover 



the o bserved positron fraction between 1 and 100 GeV. ATIC's ( jChang et al. 



20081 ) electron excess at ~ 600 GeV is explained, in this picture, as the contribu- 
tion of a few known nearby SNRs. The apparent coincident similarity between 
the cooling time of electrons at 10 GeV (where the positron/electron ratio up- 
turn), ~ 10 Myr, and the CRs protons cosmogenic age at the same energy is 
predicted by this model. 



PAMELA (lAdriani et al. 1120091 ) discovered that the CR positron/electron ratio increases 
with energy above ~ 10 GeV. This ratio should decrease according to the standard sce- 
nario, in which CR positrons are secondaries for med by interactions betwee n the primary 



CR protons and the interstellar medium (ISM) (IMoskalenko fc Strong! 1 1 9981 ). This appar- 
ent discrepancy is now commonly known as the "PAMELA anomaly". It is commonly 
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excess of CR elec trons at energies of 300 — 800 GeV. At even higher energies (1 — 4 TeV) 
HESS measures ( IH. E. S. S. Collaboration! 120081 ) a sharp decay in the electron spectrum. 
ATIC's results are usually considered as support for a dark matter interpretation for the 
PAMELA anomaly, where the observed excess corresponds to the WIMP mass. 

In the standard picture, C Rs below the knee ar e thought to originate in SNR shocks. 
This is indicated by synchrotron (IKoyama et al. I ll995h and inverse-Compton (ITanimori et al. 
19981 ) emission of high ene rgy electrons in SNRs, and the 7-ray emission, which is possibly 
from high energy protons (lAharonian et al. 1 120041 ) . Theoretical models for the CR flux de- 
scribe CR propagation in the Galaxy. CRs diffuse within the disk, and escape once they reach 
the halo height, In ~ 1 kpc, above the disk. Most CR diffusion models approximate the diffu- 
sion coefficient as D = D (E / EqY and assume that CRs are produced with a power-law spec- 
trum, Ne = dN/dE oc E~ a . The observed spectrum is then a convolution of the source spec- 
trum and propagation losses, giving for the primary electrons <fi~(E) oc E~( ae+l3 \ Positrons 
are secondary CRs formed from CR protons, and suffer additional propagation loses, implying 
(f> + (E) oc <pp(E)E~ 13 oc E~( ap+2/3 \ where ± and <p p are the CR positrons, electrons and pro- 
tons observed fluxes. The predicted flux ratio is (fi + /(</>" + (f> + ) m <p + /<fr~ oc E Qe_ar ^, where 
a e and a p are the source power-l aw indices of electrons and protons respectively. Both elec- 
trons and protons are expected (IBlandford fc Eichlerl 119871 ) to have similar spectral slopes, 



i.e., a e ~ a p , which is somewhat larger than 2. This is also supported by synchrotron rad i- 



ation observed from SNRs, which confirms the slope for the electrons ( IDuric et al. Ill995l ). 
Consequently, a p — a e < (3 ~ 0.3 — 0.6 and the standard model predicts, in contrast to 
PAMELA observations, a decreasing (f) + / <fi~ . 

The diffusing electrons and positrons cool via synchrotron and inverse-Compton scat- 
tering, with dE/dt = —bE 2 . This steepens both the electron and positron spectra at an 
energy where the cooling time equals the typical electron and positron age. However, since 
both suffer the same loses, this does not affect <p + / <p~ . Additional effects such as spallation 
and annihilation can be safely ignored at the energies of interest. 

This standard model assumes a homogenous, or at least a smoothly varying (on a 
galactic scale), source distribution (IMoskalenko fc Strongill998l ; IStrong fc Moskalenkdll9 98'l. 
However, since in s piral galaxies star formation is concentrated in spiral arms (ILacey fc Duric 
200ll ; IShavivll2003l ) and SNRs are the canonical sources of CRs, one should consider the effect 
of inhomogeneities in the CR source distribution on intermediate scales (i.e., scales smaller 
than the Galactic size but large enough such that discrete sources do not have a strong 
effect) on the CR spectrum. This inhomogeneity of sources influences the electrons/positrons 
spectra via cooling which sets a typical distance scale that an electron/positron with a given 
energy can diffuse away from its source. For a homogenous distribution cooling affects 
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the spectra of (primary) electrons and (secondary) positrons in the same way and their 
ratio is unaffected. On the other hand, primary electrons will be strongly affected by an 
inhomogeneous source distribution at energies for which the diffusion time is longer than 
the cooling time. Protons are not affected by cooling and are therefore distributed rather 
smoothly in the galaxy even if their sources are inhomogeneous. The secondary positrons 
(that are produced by the smoothly distributed protons) are only weakly affected by the 
inhomogeneity of the sources. This effect would induce an observed signature on 
with similar properties to the one observed by PAMELA. 

Motivated by this expectation we construct, first, a simple analytic model for diffusion 
from an inhomogeneous source. Consider a source at a distance d from Earth. We model the 
solar neighborhood of the galaxy as a two dimensional slab (see fig. f). The x coordinate 
(the Galactic plane) is infinite and the y coordinate (the disk height) is finite, Ih- The source 
is at the origin and Earth is at (d, 0). A CR diffuses within this slab with a constant diffusion 
coefficient D(E), and it escapes once \y\ > Ih- The contribution of CR protons that were 
generated at time if to the flux at time to can be approximated as0: 

<f> p (d,t') oc -Lexp[-(t/r e ) - (r d /2t)}, (f) 



where t = to — t',T e ~l H /D is the typical escape time and r d rs d 2 /D is the typical diffusion 
time from the source to Earth. Integration over t for a steady source, yields: 

<f) p (d) oc -=^xp -^2r d /T e , (2) 

with a similar energy dependence (via D) as for uniformly distributed sources. The average 
age of an observed proton is a = IhQ-h + V%d)/2D m max{r e , {r e Td) x ' 2 }. 

We approximate the cooling effect on the electron's flux as <p~(d, if) oc <f) p (d, tf) exp[— t/r c ), 
where r c is the typical cooling time. Integration over t reads: 



exp 



-2\/Td/T c + T d / 



(f)-(d) oc L - ± (3) 



Dy/1+T e /T t 

If t c < minjrrf, (t"^) 1 / 2 } the electron flux drops exponentially with decreasing r c , while for 
larger r c the electron flux is proportional to D~ l (relative to the source's spectrum). This 
is different than the case of uniformly distributed sources, which shows a shallower break at 



1 We assume for simplicity that the diffusion is one dimensional. This results with an exponent once 
integrated. Two dimensional diffusion (from a linear spiral arm) would give a less transparent Bessel function. 



-4- 



r c w r e from (p (r c > r e ) oc r e oc D 1 oc E 13 into (r c < r e ) oc r c oc E 1 , both relative to 
the source's spectrum. 

The positron source function is approximately proportional to <fi p (d). As positrons and 
electrons have the same cooling rate, a source at x' contributes to the positron flux at d 
approximately <p~(x' — d). Therefore: 



t. 



4> + {d) oc I (p p (x')(p (V — d)dx ! oc ( exp 





l2r d 


exp 










For r c ^> r e , the energy dependence of + relative to the source spectrum, cpf, is 
D~ 2 oc E~ 213 while for r c <C r e , oc r c /2) oc E^ 13 ^ 1 . This behavior is similar to the one 

from uniformly distributed sources. 

Eqs. |3] and H] show that for a source at a distance d from Earth, a turnover in is 
observed at E b which satisfies r c (E b ) pa mm{r x (E b ), (T e (Eb)r x (Ei ) )) 1 / 2 } . <p + / for E < E b 
decreases, while it increases for E > E b . At the same time the typical age of CR protons 
with energy E b is a ~ max{r e , (Tei^) 1 / 2 }. Therefore a natural prediction of the model is 
a(E b ) > r c (E b ) and a comparison of the two observables can be used as a consistency test 
for the model. Moreover, over a wide range of the parameter space for which d > Ih, the 
model predicts a(E b ) ~ r c (E b ) regardless of the value of the diffusion coefficient D. 



Electrons and positrons in the ISM cool as dE/dt = —bE 2 where ( iKobayashi et al. 



2004J) b w 1.8 X lO-^GeVV 1 at 1 GeV (and b w 1.4 X lO-^GeV - ^ -1 at 1 TeV), implying 
a cooling time r c = 1/(622) pa 17 Myr at E pa 10 GeV. Observational constraints on 
the typical proton CR age are measur ed at a few 1 00 MeV. Typical age s obtained are 
18+5 Mvr Jwiedenbeck fc Greinerl \imh , 27^ 9 Myr jLukasiak a/. I Il994f ) or 30l?J Myr 
( Simpson &: Garcia-Munozl 19881 ). At 10 GeV, the age should be smaller by a factor of 
~ 1 — 3, depending on the exact energy dependance of the diffusivity. Thus, according to the 
observations a(10GeV) ~ r c (10GeV) ~ lOMyr. This apparent coincidence which is explained 
naturally by our model encourages us to look for a dominant CR source at a distance of a 
~ kpc from earth. Indeed, the nearest spiral arm to Earth is the Sagittarius-Carina arm at 
a distance of ~ 1 kpc, which is just the distance needed to explain PAMELA'S observations. 

To demonstrate quantitatively the potential of this model to recover the observed be- 
havior of + /0~, we simulated num erically the C R diffusion for a realistic spiral-arm concen- 
trated source distribution (see also IShavivl 120031 ) . Before presenting these results we stress 
that all other models explaining PAMELA invoke a new ad hoc source of high energy CR 
positrons which has a negligible effect on low energy CR components. However, in our 
model, the PAMELA explanation is intimately related to low and intermediate energy CR 
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propagation in the Galaxy. Namely, by revising the source distribution of CRs, we affect 
numerous properties of ~ GeV CRs. Given that the interpretation of observations (in par- 
ticular, isotopic ratios) used to infer model parameters (such as D , (3 or l H ) depend on the 
complete model, one should proceed while baring in mind that these parameters may differ 
in our model from present canonical values. In this sense, the objective of this letter is not to 
carry a comprehensive parameter study, fitting the whole CR data set to an inhomogeneous 
source distribution model. Instead, our goal is to demonstrate the potential of the model to 
explain naturally the PAMELA anomaly. To this end we use the simplest possible model, 
fixing all parameters with the exception of the halo size, Ih, and the normalization of the 
diffusion coefficient, D Q , that we vary to fit the data. 

The geometry of the model is described in fig. [TJ We assume a spiral arm/disk SNe 
ratio of 10. The overall normalization of the sources was fit to give the electron spec- 
trum at 10 GeV. The positron production was normalized to give the positron to electron 
ratio at the same energy. For the ISM density we took the functional dependen ce from 
Strong fc Moskalenkd (119981 ) . More on the choice of the parameters can be found in IShaviv 
j200.4 



We take a diffusivity of the form D = D (E/1 GeV) 13 for E > 4 GeV and D = 
Dq(4 GeV/1 GeV)^ for E < 4 GeV. It was r ealized that such a break is required to ex- 
plain the observed break in the CR B/C ratio Jstrong fc Moskalenkd [l998l ^ (though it does 
not play an important role here). We take (3=1/3 (corresponding to turbulence with a 
Kolmogorov spectrum) and a e = a p = 2.37 such that the predicted proton spectrum will be 
consistent with the observed proton CR slope of 2.7. We also take D = 6 x 10 27 cm 2 /sec, 
which reproduces the break energy in the electron spectrum and the positron fraction. As 
predicted by the analytic model the cosmogenic age we obtain in the simulation (14 Myr at 1 
GeV per nucleon) is consistent with the observations, without fitting for it. Not surprisingly, 
the halo size and diffusivity considered here are somewhat different (on the low side) relative 
to standard values often found in the homogenous model. 

Small scale inhomogeneities are important at energies larger than a few hundreds GeV, 

for which the lifetime, and therefore propagation distance, of electrons is so short that the 

electron spectrum is dominated by a single, or at most a few nearby sources (lAtoyan. Aharonian &: Volk 
19951 ; iKobayashi et al. 1120041 ; |Profumoll2008l ). To take this effect into account we truncate 
the "homogeneous" disk component at r < 0.5 kpc and age less than t < 0.5 Myr, and we 
add all SNRs within this 4-volume: Geminga, Monogem, Vela, Loop I and the Cygnus Loop, 
as discrete instantaneous sources. These sources were described using the analytical solution 
( lAtoyan. Aharonian fc Volklll995l ) for the diffusion and cooling from an instantaneous point 
source. For the overall normalization of the point sources, we use the synchrotron observa- 
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Fig. 1. — The galaxy is modeled as a slab of width 2Z#, with Ijj — 1 kpc, inside of which the 
CR components diffuse. Beyond y = ±Ih, the CRs escape at a negligible time. CR sources 
are located in both cylinder shaped arms with a Gaussian cross-section of width o = 300 pc, 
and disk sources, with a vertical scale height of 100 pc. The assumption of straight cylinders 
is permissible given the small spiral arm pitch angle. This also makes the problem effectively 
two dim ensional. We model the Milky Way as having four spiral arms, with a pitch angle of 
% ~ 15° (j Vailed 120081 ) . implying that the arm separation (in the direction perpendicular to 
the arm axis) is d ~ (tt /2)R Q smi m 3 kpc, while the Sun is at a distance x ~ 1 kpc from 
the nearest spiral arm. Due to the motion of the arms, there is a small drift term carryin g 
the CRs away from them. For a spiral arm periodicity of P s ~ 150 Myr (jShavivl 120031 ). 
one obtains a velocity of v s ~ (it / 2) (R & sin i / P s ) « 20 km/s, which is slower than the two 
comparable diffusion times In/^e ~ x/t x w 100 km/s. A second component resides in the 
disk, with an exponential vertical decay. Because nearby sources are considered, the The 
smooth disk distribution is truncated for r < 0.5 kpc and t < 0.5 Myr. 
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tions of SN1006, which togethe r with the X-rays constrain the total energy and magnetic 



field ( lYoshida &: Yanagitalll997l ). In particular, electrons with energy > 1 GeV are found to 
carry » 2 x 10 48 erg, corresponding to 0.2% out of the total ~ 10 51 erg mechanical energy 
in SNRs. We assume that all nearby sources are similar. Note that due to their very young 
age, the discrete sources contribute a negligible amount of positrons, nor do they offset the 
cosmogenic age. 

The lower panel of fig. 2 depicts + /(0 + + obtained by the simulation. As expected 
from the simple analytical model, the fraction decreases up to ~ 10 GeV and then it starts 
increasing. This explains the so called PAMELA anomaly. As the CR protons and antiproton 
spectra are unaffected our results are consistent with PAM ELA'S observations of no excess 
in the anti-proton/proton ratio at the same energy range (jAdriani et al. 2009 ). At about 



100 GeV, the ratio flattens and it decreases above this energy because of the injection of 
"fresh" CRs from recent nearby SNRs whose high energy primary electrons don't have time 
to cool. These sources also contribute to higher energy electrons detected by ATIC. 

The upper panel of fig. 2 depicts the electron spectrum and its constituents — primary 
spiral arm electrons, primary disk electrons (without nearby sources), nearby sources and 
secondary pairs. There are two bumps in the E 3 Ne plot. The lower energy bump arises 
from spiral arm electrons, the higher energy of which cannot reach us due to cooling. The 
higher energy bump, which corresponds to the ATIC peak, is due to a few nearby SNRs. The 
three "steps" are due to the cooling cutoffs from Geminga, Loop I and the Monogem SNRs. 
Note that the high energy behavior is very sensitive to the exact diffusion model parameters 
and the poorly constrained SNR energy output in electrons. Thus, the ATIC peak is not a 
prediction of the model but rather one possible outcome. 

While the predictions for for the spiral arms CR model are very different than for 

a homogenous sources distribution, the effect on the electron spectrum is much more subtle. 
Both models predict a break of the electron spectrum at 10 GeV, The break predicted by 
spiral arm model is from a power law to an exponential, while in the homogenous model it is 
a broken power- law. Given that above ~ 100 GeV the electron spectrum is strongly affected 
by the sources that produces the ATIC bump (e.g. local SNRs), the energy range between 
10 to 100 GeV is too short to distinguish, based on the electron spectrum alone, between 
the two models. Thus, while both models can adequately reproduce the observed electron 
spectrum (at least up to 100 GeV), only the inhomogeneous source model can explain the 
observed 

One of the interesting predictions of the model where the ATIC peak is explained as 
consequences of propagation effects from local SNRs, is that the electron spectrum around the 
ATIC peak is dominated by nearby sources. These source produce only primary electrons and 
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Fig. 1.— Bottom Panel: Model results and the measured PAMELA points for the 
positron fraction. The shaded region is the variability expected from solar modulation effects 
( jClem et al. Ill996l ). Top Panel: The expected electron and positron spectra - Primary arm 



electrons (long dashed purple), primary disk electrons with nearby sources excluded (short 
dashed green), nearby SNRs (dot-dashed black), secondary positrons (dot-dashed red), and 
their sum (blue). The hatched region describes the s olar modulation range ( from 200 MV to 
1200 MV). The thre e data sets plotted are o f HEAT JDuVernois et ZlboOlh (circles), ATIC 
(jChang et al. 1120081 ) (triangles) and HESS (IH. E. S. S. Collaboration! 120081 ) (open squares). 
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have only negligible contribution to secondary positron flux. As a result ATIC observations 
force the electron/positron ratio to start decreasing at a few hundred GeV, which is not far 
above the present PAMELA measurement. It should reach a minimum around the ATIC 
peak, where it should start rising again. Whether or not it can go up to about 50% at a few 
TeV depends on whether the CRs from very recent SNe, the Cygnus Loop and Vela, could 
have reached us or not. This critically depends on the exact diffusion coefficient. Here it is 
also worth pointing out that above a few TeV the secondaries must be produced within the 
local bubble, implying that their normalization should be ten times lower than for the lower 
energy secondaries. These predictions are in contrast to the case where the ATIC peak is due 
to a primary source of pairs, in which case the positron fraction is expected to keep rising 
also at a few hundreds GeV. With these predictions, it will be straightforward in the future 
to distinguish between propagation induced "anomalies", and real anomalies arising from 
primary pairs (in particular, when PAMELA'S observations will extend to higher energies). 
Of course, it is possible that the ATIC peak is due to a source of primary pairs, while the 
PAMELA anomaly is a result of SNRs in the spiral arms, but then it would force us to 
abandon the simplicity of the model, that the anomalies are all due to propagation effects 
from a source distribution borne from the known structure of the Milky Way. 

Irrespectively, this work demonstrates that the intermediate scale inhomogeneities ex- 
pected in the CR source distribution leave nontrivial imprints on the electron and positron 
spectra. These should be further investigated before reaching definitive conclusions about 
the existence of primary positron sources. 

We thank Marc Kamionkowski, Re'em Sari and Vasiliki Pavlidou for helpful discussions. 
The work was partially supported by the ISF center for High Energy Astrophysics, an ISF 
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